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Some Probabilistic Aspects 
of Set Partitions 


Jim Pitman 


1. INTRODUCTION. A partition of the set N (1 == {1, 2, , n) is an unordered 
collection of non-empty subsets of N n . Let P„ denote the set of all such partitions, 
and let B n = #(P„), the number of partitions of PJ„. The numbers B n are known 
as Bell numbers after E.T. Bell [3, 4, 45]. See Rota [50] and Gardner [24, Chapter 2] 
for surveys of their properties and applications. The remarkable Dobihski formula 
[18] 


B n = e~'Z - 7 (n = 1 , 2 ,...) ( 1 ) 

m = l 

leads [36, 1.9] to the asymptotic evaluation 

B n j=r \(n) n + 1/2 e Kn) ~ n ~ l as n -> oo, (2) 

\n 

where A(n)log(A(n)) = n. As noted by Comtet [11], for each n the infinite sum in 
( 1 ) can be evaluated as the least integer greater than the sum of the first 2 n terms. 

From a probabilistic perspective, the series on the right side of Dobinski’s 
formula represents the nth moment of the Poisson distribution with mean 1. So 
the initially surprising fact that this series yields an integer for all n amounts to the 
fact that all positive integer moments of the Poisson(l) distribution are integers. As 
explained in Section 2, Dobiriski’s formula reduces to the fact that the factorial 
moments of the Poisson(l) distribution are identically equal to 1, and this identity 
can be understood probabilistically with essentially no calculation. 

While such probabilistic interpretations of identities related to set partitions are 
the main theme of this paper, Section 1.2 recalls an elementary combinatorial 
proof of Dobinski’s formula. 

1.1 Notation. Following the notation of [27], let denote the number of 
partitions of N n into exactly k distinct non-empty subsets, so that 

*-£{ 2 }- (3) 

The j" j are known as the Stirling numbers of the second kind. Let m- denote the 
falling factorial with k factors 

m- = m(m - 1) ••• (m — k + 1), (4) 

which, for positive integers m and k , is the number of permutations of length k of 
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m distinct symbols. The formula 



( 5 ) 


decomposes the number m n of sequences of m distinct symbols of length n as the 
sum over k of the number of such sequences that contain exactly k distinct 
symbols [54, p. 35]. As an identity of polynomials in m of degree n , this identity 

provides an alternative definition of the coefficients for 1 < k < n. See 

[11, 47, 48, 54] for background and a wealth of further information about Stirling 
numbers. 


1.2 A quick proof of Dobinski’s formula. This argument is attributed to Schiitzen- 
berger by Foata [22, p. 73]. Divide (5) by m! to obtain for positive integers m 
and n 


m n 

m\ 



1 

(m — k)\ 


( 6 ) 


This is the identity of coefficients of A m in the power series identity 



( 7 ) 


which, upon rearrangement, gives the following horizontal generating function for 
the Stirling numbers of the second kind: 



Now take A = 1 and use (3) to obtain Dobinski’s formula (1). The polynomial 
appearing in (8) is known as an exponential polynomial Many other proofs of the 
generalization (8) of Dobiriski’s formula are known. See for instance Roman [49, 
p. 66] and Wilf [58, p. 106]. Closely related arguments appear in Rota [50], Berge 
[5, p. 44], Comtet [11, p. 211], Lupas [37], and Chen-Yeh [10]. 


2. MOMENTS. For a non-negative integer-valued random variable X with 

P{X = tn)=p m (m = 0,1,...) (9) 

and a non-negative function /, let 

E[f(X)] ■■= (10) 

m 

which is the expected value of f(X) for X with distribution (9). See [20, 43] for 
background. From (5) and linearity of the expectation operator E , we obtain the 
following well-known formula for E[X n ], the rcth moment of X, in terms of E[X~], 
the ktYi factorial moment of X for 1 < k < n [47, 14]: 

*[*"]-£{;;}*[**]• (ii) 

For A > 0, let X x denote a random variable with the Poisson distribution 

A m 

P{X K = m) = e- A — (m = 0, 1, . . . ) (12) 

m ! 
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so that 


oo 

E[f{X A )] = <T A E/(W) — . (13) 

m = 0 m • 

Take /(m) = m" to see that the right side of (8) equals E[X"\ so the identity (8) 
amounts to the formula 

E[X n=p i { n k )* lc ( n = (14) 

for the moments of the Poisson(A) distribution [46,42]. This formula is the 
particular case of (11) for X with Poisson(A) distribution, for it is known [46, 14] 
that 

E[xf\= A* (k = 1 , 2 , ... ). (15) 

Formula (15) follows easily from (13) with f(m) = m- by change of summation 
variable from m to j = m — k. In particular, for A = 1 the factorial moments of 
the Poisson(l) distribution are identically equal to 1. So Dobinski’s formula (1) can 
be read from (14) for A = 1, which follows as indicated above from (11) and (15). 
In essence, this is Rota’s [50] proof of Dobinski’s formula cast in probabilistic 
notation. This argument differs from the proof in Section 1.2 in that it involves 
checking (15) for A = 1. 

Formula (15) has the following interpretation in terms of a Poisson process 
[33, 43]. Let 

0 < U a) < - < U m < 1 (16) 

denote the random locations in (0, 1) of the points of a homogeneous Poisson 
process on (0, 1) with mean intensity measure A du for 0 < u < 1. For each 
k = 1,2,... define an associated k-tuple point process , with points in (0, l) k , to 
have a point at each of the locations (U ((Ti) , . . . , U ((Tk) ) as a ranges over the X£ 
different permutations of {1, ... , X k ) of length k. For distinct u t e (0, 1), indepen- 
dence properties of the basic Poisson process on (0, 1) imply that the mean 
intensity of the k - tuple point process at (u l9 . . . , u k ) e (0, 1)^ is 

P(some e du t for each 1 < i < k) (A du x ) ••• (A du k ) 

du x ••• du k du 1 ••• du k ’ ^ ^ 

so the expected number of points in the A-tuple point process is 



Constantine and Savits [12] derive a generalization of Dobinski’s formula by 
consideration of compound nonhomogeneous Poisson processes. See also Stam 
[52] and Di Bucchianico [8] for related results. For various applications of Stirling 
numbers and their generalizations to the computation of moments of probability 
distributions, see [47, 9]. Moments of the normal distribution also have interesting 
combinatorial interpretations [19,25]. More generally, the idea of representing 
combinatorially defined numbers by an infinite sum or an integral, typically with a 
probabilistic interpretation, has proved to be a very fruitful one. Other examples 
are the representation of n\ as a gamma integral, which leads to Stirling’s formula 
[7,16,38], and Laplace’s representation of £:th differences of powers [35,14,30], 
which yields an asymptotic formula for the Stirling numbers of the second kind. 
See [41] for a recent survey of asymptotic enumeration methods. 
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3. VARIATIONS OF DOBINSKI’S FORMULA. The derivation of Dobinski’s for- 
mula given in the previous section yields the following proposition: 


Proposition 1. Let X be a random variable with values in {0, 1, 2, ... } and let n be a 
positive integer. The following two conditions are equivalent : 

(i) the first n factorial moments of X are identically equal to 1; 

(ii) the kth moment of X equals B k for every 1 < k < n. 

It is well known that for each A > 0 the Poisson(A) distribution is uniquely 
determined by its moments; see for instance [6, Section 30]. The Poisson(l) 
distribution is therefore the unique probability distribution whose nth moment 
equals B n for every n. But for each fixed n there are many probability distribu- 
tions on {0, 1, 2, ... } that have the same first n moments as Poisson(l). It is obvious 
that there can be at most one such distribution of X with P(X <n) = 1, because 
the moment conditions amount to a system of n linearly independent equations in 
n unknowns Less obvious is the fact that the unique solution of these 

equations is such that p t > 0 for 1 < i < n and EJL x p t < 1, so that (p l9 . . . , p n ) is 
the restriction to {1, ... , n } of a unique probability distribution on {0, 1, ... , n). But 
this probability distribution on {0, 1 , 2, ... , n }, whose first n factorial moments are 
identically equal to one, is known to arise in the setting of the classical matching 
problem [31,14,20,56]. If M n is the number of fixed points of a uniformly 
distributed random permutation of N n , then it is easy to show by the method of 
indicators that the first n factorial moments of M n are identically equal to 1; see 
[14]. The distribution of a random variable X with range {0, 1, ... , n) is recovered 
from its factorial moments by the classical sieve formula [14] 

1 * { — l) mk E\X-\ 

P(X = m) = - L ~ T 777 — “ (m = 0,l, (19) 

m\ k=m ( rn-k)\ 

For X = M n with E[M%] = 1 for 0 < k < n, this simplifies to 

1 (- 1 )* 

P(M n = m) = — - £ j — (m = 0,1,..., ra). (20) 

ml s=0 s\ 


See [20, Section IV.4] for further discussion. According to Proposition 1, the kth 
moment of M n equals B k for every 1 < k < n. That is to say, 


B k = E 


m 


k n—n 


(~iy 


m * 

m = \ m * j— 0 


S ! 


(1 < k < n). 


( 21 ) 


This variation of Dobinski’s formula is derived in quite a different way by Wilf [58, 
p. 22] by substituting the classical formula 



1 X 

Km 7 = 1 


k~j 


( 22 ) 


into (3). As observed by Wilf, Dobinski’s formula (1) follows easily from (21) by 
letting n -> oo. See also Lovasz [36, 1.9] for a similar argument, James and Kerber 
[32, pp. 227-237] for connections with the representation theory of the symmetric 
group, and Diaconis and Shashahani [17] for various generalizations. In Dale and 
Skau [13] the Bell numbers appear as the factorial moments of a probability 
distribution on the non-negative integers. 
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4. THE MOMENT GENERATING FUNCTION. Consider now the 


ating function ( m.g.f .) of the Poisson(A) distribution: 


E[exp(0X A )] =E 


E 

n = 0 


0 n Xf 


n\ 


E £[*;] 


n = 0 


0 " 



moment gener- 


(23) 


where the series converge for all real 6 and the interchange of E and E is easily 
justified. See [6] for a modern treatment of m.g.f s. From (13) with /(m) = e 6m 
there is the standard formula 


( a v rn 

Xe ® ) 

£.[exp(WA A jj =e ■' L — = exp(A(e e - 1)). (24) 

m = 0 m ‘ 

This combines with (8) to yield the following double generating function of the 
Stirling numbers of the second kind. This classical formula [11, p. 206] is an 
identity between two different expressions for the m.g.f. in 6 of the Poisson(A) 
distribution: 

00 n ( \ x k d n 

1 + „? 1 ,? 1 {^}'^r = exp(A(e9 " 1)) - (25) 

In particular, for A = 1 this reduces by (3) to Bell’s [3, 4] formula 

oo 0 n 

1 + E B n — = exp(e 5 * * * 9 - 1), (26) 

which gives two expressions for the m.g.f. in 0 of the Poisson(l) distribution. 
Equating coefficients of X k in (25) yields the vertical generating function of the 
Stirling numbers of the second kind: 

£{:}£= < 2T > 

See [11,47,54] for alternative derivations of these identities. There are similar 
identities for many other arrays of combinatorial numbers, such as the binomial 
coefficients and Stirling numbers of the first kind [11,58], [27, p. 351], most of 
which admit probabilistic interpretations. Formulae with binomial coefficients 
typically involve independent trials, while those with Stirling numbers of the first 
kind typically involve the cycle structure of random permutations [1]. See also [2] 
for probabilistic analysis of more general combinatorial structures and further 
references. 


5. RANDOM PARTITIONS. A random partition of N n is a random variable II 
with values in the set P n of partitions of N n . The distribution of II then refers to 
the collection of probabilities P(Il = it) as tt ranges over P n . Questions about 

enumeration of partitions of of various kinds can be phrased probabilistically in 

terms of a uniform random partition , that is, a random partition II with the 
uniform distribution P(U = tt) = 1 /B n for each partition tt P n . For develop- 
ments of this idea see [29,28,51,23]. Random partitions with non-uniform distri- 
bution also arise naturally in various contexts, so it is useful to have models for 
random partitions, both uniform and non-uniform. 

The following random allocation scheme provides a basic method of generating 

a random partition of See [14,34,57] for extensive study of this and related 
schemes, and further references. Throw n balls labelled by into m boxes 
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labelled by f\l m , and assume that all m n possible allocations of balls into boxes are 
equally likely. Let 1 1 „ be the partition of balls by boxes. More formally, let X t be 
the number of the box containing the ith ball for 1 < / < n. Then the are 
independent and uniformly distributed on N ffl , and II „ m is the partition of 
induced by the random equivalence relation i ~ j if and only if X t = X,. Formally, 
the Xj can be regarded as coordinate maps defined on (N„) m , and II„ m is then 
defined as a map from (N„) m to P„, the set of partitions of N„. Let #( 7 r) denote 
the number of subsets in a partition ire P n . The distribution of 1I„„, induced by 
the uniform distribution P on can be read from formula (5): 

m- 

p( n„ m = tt) = — if #(tt) = k. (28) 

The distribution of #(II„ m ), the number of occupied boxes when n balls are 
thrown into m boxes, is given by the following probabilistic equivalent of (5): 

P[#(U nm )=k] = { n k }^ (1 < k < n). (29) 

Because the probability displayed in (28) depends on the number of occupied 
boxes k, for n > 3 this random partition II of does not have uniform 
distribution on P n for any m. However, as observed by Stam [53], for each fixed n 
it is possible to generate a uniformly distributed random element of P n by a 
suitable randomization of m. The following proposition was suggested by Stam’s 
construction, which is described in Corollary 3. 

Proposition 2. Let M be a random variable with values in {1, 2, . . . }, and suppose 
given M = m that n balls labelled by N n are thrown independently and uniformly at 
random into m boxes. Let U nM denote the random partition of N n so generated. The 
following two conditions are equivalent'. 

(i) U nM has uniform distribution over the set P n of all partitions of N n ; 

(ii) the distribution of M is of the form 

m n p m 

P(M = m) = — " (m = 1,2,...) (30) 

B n 

for some probability distribution (p m ) on {0,1,2,...} whose first n factorial 
moments are identically equal to 1. 

Before the proof, here are two corollaries which follow immediately from the 
Proposition and the discussion in Sections 2 and 3: 

Corollary 3. [53] If M has the distribution (30) for p m = e~ x /m\, then U nM has 
uniform distribution on 

Corollary 4. For each n there is a unique distribution of M such that 

P(M < n) = 1 and U nM has uniform distribution on P n . 

This distribution of M is defined by (30) for p m = P(M n = m) as in (20) with M n the 
number of fixed points of a uniform random permutation of N n . 
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( 31 ) 


Proof of Proposition 2. By conditioning on M and using (28), 

“ m- 

P(U nM = 77 ) = £ — P(M = m) if #(7T) = fe, 

m = 1 m 

so the distribution of TV nM is uniform on ?,. if and only if 

" m- 1 

E-^(M = m) = — (1 Zk<n). (32) 

m = 1 m B n 

Define 

Pm =B n m~ n P(M = m) (m = 1,2,...), (33) 

so that (32) becomes 

oo 

£ = 1 (1 <k<n), (34) 

m = 1 

which for k = 1 implies that £“ =1 p m < YT m=l mp m = 1. It follows that U nM is 
uniform if and only if ( p m ) derived from the distribution of M via (33) is the 
restriction to (1, 2, ... } of a probability distribution on (0, 1, 2, ... } whose first n 
factorial moments are equal to 1. This is condition (ii). ■ 

As shown by Stam, Corollary 3 allows numerous results regarding the asymp- 
totic distribution for large n of a uniform random partition of to be deduced 
from corresponding results for the classical occupancy problem defined by random 
allocations of balls in boxes, for which see [34, 57]. See also [2, 15, 23, 26, 28, 29, 51] 
for a more detailed account of the asymptotics of uniform random partitions of 
As a variation, the following corollary is easily obtained by a similar argument: 

Corollary 5. Suppose that M has the distribution 

m n P(X x = m) 

P(M = m) = (m = 1,2,...), (35) 

^n{ A ) 

where X K has the PoissonOO distribution (12), and /x,,( A) = E(X" ). Then the 
distribution of II „ M is given by 

A* 

P(U nM = it) = — — - if#(n)=k. (36) 

Pn\ 

As a check, (36) implies 

p[#(n nM ) = k] = (1 <k<n). (37) 

The fact that these probabilities sum to 1 amounts to formula (14) for p n ( A). The 
distribution of II„ M defined by formula (36) defines a Gibbs distribution on 
partitions of N n . See [55,44] for further discussion of such Gibbs distributions on 
sets of combinatorial objects. See Nijenhuis and Wilf [39] for a recursive algorithm 
to construct a uniform random partition of based on the recurrence 

B n -1 + (38) 

where the right side counts the number of partitions it of N n according to the size 
k of the subset in tt that contains n [36, Problem 1.10]. See [40] for related 
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combinatorial algorithms, and [21] for a recent systematic approach to the random 
generation of labelled combinatorial structures and further references on this 
topic. 


ACKNOWLEDGMENTS. Thanks to Persi Diaconis, Hsien-Kuei Hwang, and the referee for providing 
some references. 


REFERENCES 


1. R. Arratia and S. Tavare. The cycle structure of random permutations, Ann. Prob., 20:1567-1591, 
1992. 

2. R. Arratia and S. Tavare. Independent process approximations for random combinatorial struc- 
tures, Adv. Math., 104:90-154, 1994. 

3. E. T. Bell. Exponential numbers, Amer. Math. Monthly , 41:411-419, 1934. 

4. E. T. Bell. Exponential polynomials, Annals of Mathematics, 35:258-277, 1934. 

5. C. Berge. Principles of Combinatorics. Academic Press, New York and London, 1971. 

6. P. Billingsley. Probability and Measure, 2nd ed. Wiley, New York, 1986. 

7. C. R. Blyth and P. K. Pathak. A note on easy proofs of Stirling’s theorem, Amer. Math. Monthly, 
93:376-379, 1986. 

8. A. Di Bucchianico. Representations of Sheffer polynomials, Studies in Applied Math., 93:1-14, 
1994. 

9. Ch. A. Charalambides and J. Singh. A review of the Stirling numbers, their generalizations and 
statistical applications, Commun. Statist. -Theory Meth., 17:2533-2595, 1988. 

10. B. Chen and Y.-N. Yeh. Some explanations of Dobinski’s formula, Studies in Applied Math., 
92:191-199, 1994. 

11. L. Comtet. Advanced Combinatorics. D. Reidel Pub. Co., Boston, 1974. (translated from French). 

12. G. M. Constantine and T. H. Savits. A stochastic representation of partition identities, SIAM J. 
Discrete Math., 7:194-202, 1994. 

13. K. Dale and I. Skau. The (generalized) secretary’s packet problem and the Bell numbers, Discrete 
Math., 137:357-360, 1995. 

14. F. N. David and D. E. Barton. Combinatorial Chance. Griffins, London, 1962. 

15. J. M. DeLaurentis and B. G. Pittel. Counting subsets of the random partition and the ‘Brownian 
bridge’ process, Stochastic Processes and their Applications, 15:155-167, 1983. 

16. P. Diaconis and D. Freedman. An elementary proof of Stirling’s formula, Amer. Math. Monthly, 
93:123-125, 1986. 

17. P. Diaconis and M. Shahshahani. On the eigenvalues of random matrices, J. Applied Probability, 
Special vol. 31A:49-62, 1994. 

18. G. Dobinski. Summirung der Reihe Ln m /n\ fiir m = 1, 2, 3, 4, 5, . . . , Grunert Archiv {Arch, fur 
Mat. und Physik), 61:333-336, 1877. 

19. E. B. Dynkin. Gaussian and non-Gaussian random fields associated with Markov processes, 
J. Funct. Anal., 55:344-376, 1984. 

20. W. Feller. An Introduction to Probability Theory and its Applications, Vol. 1, 3rd ed. Wiley, New 
York, 1968. 

21. Ph. Flajolet, P. Zimmerman, and B. V. Cutsem. A calculus for the random generation of labelled 
combinatorial structures, Theoretical Computer Science , 132:1-35, 1994. 

22. D. Foata. La serie generatrice exponentielle dans les problemes d’ enumeration, volume 54 of 
Seminaire de Mathematiques superieures. Presses de PUniversite de Montreal, 1974. 

23. B. Fristedt. The structure of partitions of large sets, Technical Report 86-154, University of 
Minnesota Mathematics Dept., Minneapolis, Minnesota, 1987. 

24. M. Gardner. Fractal music , HyperCards and more — : mathematical recreations from Scientific 
American magazine. W.H. Freeman, New York, 1992. 

25. C. D. Godsil. Algebraic Combinatorics. Chapman & Hall, New York, 1993. 

26. W. M. Y. Goh and E. Schmutz. Random set partitions, SIAM J. Discrete Math., 7:419-436, 1994. 

27. R. L. Graham, D. E. Knuth, and O. Patashnik. Concrete Mathematics: a foundation for computer 
science, 2nd ed. Addison-Wesley, Reading, Mass., 1989. 

28. J. Haigh. Random equivalence relations, Journal of Combinatorial Theory, 13:287-295, 1972. 

29. L. H. Harper. Stirling behavior is asymptotically normal, Ann. Math. Stat., 38:410-414, 1966. 

30. L. Holst. On numbers related to partitions of unlike objects and occupancy problems, Europ. J. 
Combinatorics, 2:231-237, 1981. 


208 


SOME PROBABILISTIC ASPECTS OF SET PARTITIONS 


[March 



31. J. O. Irwin. A unified derivation of some well-known frequency distribution of interest in biometry 
and statistics, J.R. Statist. Soc ., 118:389-404, 1955. 

32. G. D. James and A. Kerber. The representation theory of the symmetric group , volume 16 of 
Encyclopedia of mathematics and its applications. Addison-Wesley, Reading, Mass., 1981. 

33. J. F. C. Kingman. Poisson Processes. Clarendon Press, Oxford, 1993. 

34. V. F. Kolchin, B. A. Sevastyanov, and V. P. Christyakov. Random Allocations. V. H. Winston & 
Sons, Washington, D. C., 1978. 

35. P. S. Laplace. Theorie Analytique des Probabilities, 2nd ed.. Paris, 1814. 

36. L. Lovasz. Combinatorial Problems and Exercises , 2nd ed. North-Holland, Amsterdam, 1993. 

37. A. Lupas. Dobinski-type formula for binomial polynomials, Stud. Univ. Babes-Bolyai Math., 
33(2):40-44, 1988. 

38. V. Namias. A simple derivation of Stirling’s asymptotic series, Amer. Math. Monthly, 93:25-29, 
1986. 

39. A. Nijenhuis and H. S. Wilf. A method and two algorithms on the theory of partitions, J. Combin. 
Theory A, 18:219-222, 1975. 

40. A. Nijenhuis and H. S. Wilf. Combinatorial algorithms for computers and calculators, 2nd ed. 
Academic Press, New York, 1978. 

41. A. M. Odlyzko. Asymptotic enumeration methods, in R. L. Graham, M. Grotschel, and L. Lovasz, 
editors, Handbook of Combinatorics Vol. II, pp. 1063-1229. Elsevier, New York, 1995. 

42. C. Philipson. A note on the moments of a Poisson probability distribution, Skand. Aktuarietidskrift , 
46:243-4, 1963. 

43. J. Pitman. Probability. Springer-Verlag, New York, 1993. 

44. J. Pitman. Probabilistic bounds on the coefficients of polynomials with only real zeros, Technical 
Report 453, Dept. Statistics, U.C. Berkeley, 1996. To appear in J. Comb. Theory A. 

45. C. Reid. The search for E. T. Bell : also known as John Taine. MAA Spectrum. Mathematical 
Association of America, Washington, DC, 1993. 

46. J. Riordan. Moment recurrence relations for the binomial Poisson and hypergeometric probability 
distributions, Ann. Math. Stat., 8:103-111, 1937. 

47. J. Riordan. An Introduction to Combinatorial Analysis. Wiley, New York, 1958. 

48. J. Riordan. Combinatorial Identities. Wiley, New York, 1968. 

49. S. Roman. The Umbral Calculus. Academic Press, New York, 1984. 

50. G.-C. Rota. The number of partitions of a set, Amer. Math. Monthly, 71:498-504, 1964. 

51. V. N. Sachkov. Random partitions of sets, Theory Probab. Appl., 19:184-190, 1973. 

52. A. J. Stam. Polynomials of binomial type and compound Poisson processes, Journal of Mathemati- 
cal Analysis and Applications, 130:151, 1988. 

53. A. J. Stam. Generation of a random partition of a set by an urn model, J. Combin. Theory A, 
35:231-240, 1983. 

54. R. P. Stanley. Enumeratiue Combinatorics, Vol I. Wadsworth & Brooks/Cole, Monterey, Califor- 
nia, 1986. 

55. J. M. Steele. Gibbs’ measures on combinatorial objects and the central limit theorem for an 
exponential family of random trees, Probability in the Engineering and Informational Sciences, 
1:47-59, 1987. 

56. L. Takacs. The problem of coincidences, Arch. Hist. Exact. Sci ., 21:229-244, 1980. 

57. V. A. Vatutin and V. G. Mikhailov. Limit theorems for the number of empty cells in an 
equiprobable scheme for group allocation of particles, Theory of Probability and its Applications 
(Trans l . of Teorija Verojatnostei i ee Primenenija ), 27:734-743, 1982. 

58. H. Wilf. Generatingfunctionology. Academic Press, San Diego, 1990. 

Department of Statistics 
University of California 
367 Evans Hall # 3860 
Berkeley, CA 94720-3860 
pitman @stat.Berkeley.ED U 


1997 ] 


SOME PROBABILISTIC ASPECTS OF SET PARTITIONS 


209 



